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The flavor changing decay Z — * didj is investigated with special emphasis on the 6s final state. 
Various models for flavor violation are considered: two Higgs doublet models (2HDM's), super- 
symmetry (SUSY) with flavor violation in the up and down-type squark mass matrices and SUSY 
with flavor violation mediated by i?-parity-violating interaction. We find that, within the SUSY 
scenarios for flavor violation, the branching ratio for the decay Z — > bs can reach 10~ 6 for large 
tan/3 values, while the typical size for this branching ratio in the 2HDM's considered is about two 
orders of magnitudes smaller at best. Thus, flavor changing SUSY signatures in radiative Z decays 
such as Z — * bs may be accessible to future "Z factories" such as a Giga-Z version of the TESLA 
design. 

PACS numbers: 13.38.-b, 12.60.-i, 14.70.-e, 12.60.Jv 



I. INTRODUCTION 

Rare processes involving various particles have always 
been a gold mine for extracting interesting physics jl], ^| . 
For example, the smallness of Flavor Changing Neutral 
Currents (FCNC) in the K system prompted the intro- 
duction of the GIM mechanism and subsequently to the 
prediction m c w 1.5 GeV and the discovery of J/\f r and 
D's. BB mixing was a precursor to a heavy top quark, 
as confirmed by experiment. FCNC rare top decays, for 
which there are only weak upper bounds, will hopefully 
be discovered in future experiments, thus serving as di- 
rect indications for deviations from the Standard Model 
(SM) , since the latter leads to branching ratios which are 
smaller than 1CP 10 . 

The situation in rare Z decays, which is the subject of 
this paper, bears some similarities to rare t decays. In 
both cases the SM results from the loop induced FCNC 
decays are very small, beyond reach, at least for t, in the 
foreseeable future. Therefore, any significant detection 
of a rare decay at the level higher than 1CU 10 or 1CP 8 
for t or Z, respectively, would serve as an indisputable 
proof for physics beyond the the SM. If new physics is 
"around the corner", i.e. at « 1 TeV, the Z boson and 
the top being so close to that scale, are expected to be 
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the particles most affected by new physics. 

In this paper we study the rare decays Z — > djdj, 
where /, J indicate the generation index of a charge — 1/3 
quark, in various models. In the SM it was found that || 
Br(Z — + bs + bs) ~ 10~ 8 ; we do not repeat the SM calcu- 
lation here. Three of the models we discuss, have already 
been considered in connection with Z — > bs, namely the 
2 Higgs Doublets Model type II (2HDMII) f§, Super- 
symmetry (SUSY) § and SUSY with R-Parity Violation 
(RPV) §. Therefore we comment , wherever it is rele- 
vant in the coming sections, about differences and simi- 
larities with previous works. Note that in addition one 
can find in the literature discussion of FCNC hadronic Z 
decays, in models not covered by us in the present article 

0- 

Experimentally, the attention devoted to FCNC in Z 
decays at LEP and SLD has been close to nil. The best 
upper limit is || 

Br(Z -> bq) < 1.8 x 1(T 3 at 90% CL . (1) 

q—d.s 

This preliminary result is based on about 3.5 x 10 6 
hadronic decays, and we used |J Br(Z — > hadrons) = 0.7. 
We urge our experimental colleagues to sift through their 
LEP data to improve the current, rather loose, limit. 

In the future, there will be at least two venues in which 
Z bosons will be produced in much larger quantities than 
their number in LEP. In a high luminosity LHC with an 
integrated luminosity of 100 fb _1 , one expects 5.5 x 10 9 
Z bosons to be produced pi. A cleaner environment 
for the processes at hand, will be provided by a future 
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e + e~ linear collider. In particular, there is a viable pos- 
sibility to lower the TESLA center of mass energy down 
to y/s = mz, the so called "Giga-Z" option. With in- 
tegrated luminosity of 30 fb~ , it is possible to produce 
more than 10 9 Z bosons O], about 2 orders of magni- 
tudes larger than in LEP. To grasp the improvement in 
going from LEP to Giga-Z option of TESLA, we note 
that while the sensitivity of LEP to Z — > t/i was w 10~ 5 
0, it is expected to be ~ 10" 
Beyond LHC and the e +< 



in Giga-Z TESLA |H 
linear collider 



~e~ linear collider, there is also 
considerable interest in the community for a high energy 
muon collider |l2| . If this ever becomes a reality, it would 
also afford another very good opportunity for studying 
rare flavor changing decays and interactions [|l3| 

The paper is organized as follows: In Section 2 we 
present a generic calculation of the Zdjdj vertex at the 
one loop level. This result will assist us to evaluate the 
branching ratio for the FCNC Z decays in any partic- 
ular model. In Section 3, the results of two variants of 
the Two Higgs Doublets Model (2HDM), namely the so 
called type II 2HDM (2HDMII) and the Top-Higgs Two 
Doublets Model (T2HDM), will be reported. In 2HD- 
MII, T2HDM we get the disappointing results Br(Z -> 
bs + bs) ~ 10 -10 , 10 -8 , respectively. We then move on to 
Section 4, where our results in Supersymmetry (SUSY) 
with squark mixing, are displayed. Again, two options 
are presented, the first one with b— s mixing and the sec- 
ond one with t — c mixing. In the first case the branching 
ratio can reach a respectable Bi(Z — > bs + bs) ~ 10~ 6 
while the second case yields a branching ratio of 0(1O -8 ). 
In Section 5 we turn to SUSY with R-parity violation 
(RPV), where the effects of A' trilinear coupling terms in 
the RPV superpotential and of b terms (b is the coeffi- 
cient of the soft breaking RPV bilinear term), are con- 
sidered. Two categories of RPV are considered: Those 
which lead to a branching ratio a A' x A' and those with 
a branching ratio cx bX'. For the first category we get 
typically branching ratios at the level of 10~ 10 , while for 
the second type of RPV, we find an encouraging possi- 
bility of Br(Z -> bs + bs) ~ 10~ 6 . Finally, in Section 6 
we summarize our results. 



II. GENERIC SCALAR CALCULATION 

In this section we outline the generic framework for 
calculating the radiative one-loop flavor changing inter- 
action vertex Vdjdj with I ^ J and V = Z or 7. We 
define the one-loop amplitude for V — > didj in terms 
of form factors which are calculated for the complete 
set of one-loop diagrams that can potentially contribute 
to V — ► didj in the presence of flavor violating inter- 
actions between neutral scalars and fermions as well as 
non-diagonal vertices of charged scalars with fermions of 
different generations. 

The diagrams that modify (at one loop) the Vdjdj 
coupling due to charged or neutral scalar exchanges are 
depicted in Fig. |l]. 
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FIG. 1: One-loop diagrams that contribute to the flavor 
changing transition V — ► djdj, due to scalar- fermion ex- 
changes. 



In what follows, we will denote the internal scalars (S) 
in the loops by Greek letters and fermions (/) will be 
given Latin indices 

From Fig. [l] it is evident that we have only three types 
of interaction vertices to consider. These are defined as 
follows: 

1- Vu ~ fi ~ fj interaction 




L{Vf) J 



R 



(2) 



where L(R) = (l-(+)7 5 )/2. For the case of the SM 
couplings of a vector boson V to a pair of fermions, 
i.e., / = u (up-quark) or / = d (down-quark), there 



are only diagonal Vff couplings, 
have: 



In this case we 
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MVS)^ ~j) = a L,R(Vf) 



where 



(3) 



ticles with incoming 4-momenta p a and pp , respec- 
tively. 



3. S a — fi — dj interaction 



a-LMZf) = -9z [Tl U J - s'wQ f 
a L{if) = a R(if) = -9-yQf , 



(4) 
(5) 



with Tl (u) = 1/2 and Tp ' = -1/2 for an up and 

down-quark, respectively, and T R ' = 0. Also, Qf 
is the charge of / and 



9z 



Sy/Cw 



.97 = e 



(G) 



2. V M — S a — Sp interaction 

Sp/ 
If- 

\ 

where S a and Sp are charged or neutral spin par- 



(7) 



'a 




L(a) -R(q) 



(8) 



where d is a down-quark. 



The one-loop amplitudes M. k (k = 1,2,3,4 corre- 
sponding to diagrams 1,2,3,4 in Fig. |[ respectively), 
for the decay V — > djdj can be parametrized generically 
as follows: 



Ml J = 



16tt 



2 el{q)u dl {pi) {r (A^ k L + A% k R) + (fl£f fe L + B ! R J k R) p»} v dj (pj) , 



(9) 



where e v (q) is the polarization vector of V, q is its 4- 
momenta and (v ( j i ) is the Dirac spinor of the outgo- 
ing di with 4-momenta pi (dj with 4-momenta pj) such 
that q = pi+pj. Also, A{ J k , A R J k , B l L \, B% h are 
momentum dependent form factors. 

Using the Feynman rules in (||) , (Q) and (||) , these form 
factors can be readily calculated for each diagram. Ne- 
glecting terms of 0(rrib/y/<p) we get: 



AlJ 
A L.2 



B L,2 



- 2 E b H«AU {4(v f) [Pi 'PJ (C23 " C\ 2 ) 

-C| 4 ] • I'.l'j'iij '"./ "'/' Vj Cl) , (12) 
2 E b R(aAU {P^fAv f) (C?i ~ C\ 2 ) 



(13) 



B LA 



2 J2 P^ugv !3 b% a) b l i {p) (Cl + C\ x ) 



where Pi — — 1 if the internal fermion in the loop is a 

(10) 

charged conjugate state (/f) or else = 1. 

Combining the contribution of the two self energy dia- 
(11) grams, i.e., M3 + M4 = M34, and similarly for the form 
factors, e.g., A^ J 3 + A£ J 4 = v4£ J 34 etc., we have: 
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grams k — 1 — 4 and are given by: 



A u 



a-L(Vd) E ki(a) 6 i(a)-"l 



L,34 



= 



(14) 
(15) 



The right-handed form factors, A R 1 k and B Rk , are ob- 
tained from the corresponding left handed ones, A ! L J k and 
Bjf k respectively, by interchanging L — > i? and i? — > L 
in all the couplings which appear in (|Io|)-(15). 



The three-point one-loop form factors with x 6 
0,11,12,21,22,23,24, and two-point form factors B k x 
with x € 0,1, correspond to the loop integrals of dia- 



Cl = C x (m 2 ft ,m 2 Sa ,m 2 Sfj ,m 2 dl ,q 2 ,m 2 dj ^ , (16) 
C 2 X = C x (m 2 Sa ,m 2 ft ,m 2 fj ,m 2 dj ,q 2 ,m 2 d ^ , (17) 
Bl = B x (ml,m 2 Sa ,m dl ) , (18) 

where B x {m\,m\,p 2 ) and C x (m 2 ,m 2 ,,m 2 ,p 2 ,p 2 ,,p 2 ) are 
defined in the appendix. 

In terms of the above form factors, the partial width for 
the decay Z — > d/dj is: 



3 



16ti- 



16tt 



2 (| A? 



^ I 2 ) 



Bl I 2 



(19) 



where = 3 is the color factor and 



Al 
Bl 



AlJ 

EE B IJ 



4& 



AlJ 



PA 



B 



u 

P. 2 



B 



IJ 

P. 34 



(20) 
(21) 



for P = L and R. 



III. TWO HIGGS DOUBLETS MODELS 

In a 2HDM with flavor diagonal couplings of the neu- 
tral Higgs to down-quarks, the flavor changing decay 
Z — > djdj proceeds through the one-loop diagrams in 
Fig. | 

The interaction vertices required for the calculation of 
the form factors defined in (0) in such models are: 



V^SaSp 
Safidj 



Z^H+R- 



(22) 



H+u,d 



where the Z^UiUj coupling is the SM one as given in 
& = B + is the only charged Higgs present in 
this type of models and = itj, i = 1, 2, 3 for the three 
up-type quarks u\ = u, u 2 = c, ua = t. 

The Z^H + H~ coupling is obtained from the scalar ki- 
netic term where <I>i i2 are the two SU(2) 
Higgs doublets. This coupling is, therefore, generic to any 
version of a 2HDM. 

The coupling H + Uidj is obtained from the Yukawa po- 
tential. The most general Yukawa interaction of a 2HDM 
can be written as (see e.g., H): 



hi 



(D^ 1+ D?-$ 2 )4 



(23) 



where Ql is the SU(2) left-handed quark doublet, ur 
and dft are the right-handed up and down SU(2) singlets, 
respectively, and #i, 2 = ia 2 <£* l2 . Also, U l ,U 2 \ D l , D 2 
are general Yukawa 3x3 matrices in flavor space. The 
different types of 2HDM's are then categorized according 
to the choice of the Yukawa matrices U 1 , U 2 , D , D 2 . 

In what follows we will focus on two specific versions 
of a 2HDM: 

• 2HDM of type II (2HDMII) 

The 2HDMII follows from the choice U 1 = and 
D 2 — in which case only <E> 2 generates the masses 
of the up-type quarks and only $i is responsible for 
the mass generation of the down- type quarks [ p5[ . 
This version of a 2HDM is realized in the Minimal 
Supersymmetric Standard Model (MSSM). 

• 2HDM "for the top-quark" (T2HDM) 

In the T2HDM , the large mass of the top-quark 
is accommodated in a natural fashion by coupling 
the second Higgs doublet ($2), which has a much 
larger vacuum expectation value (VEV), only to 
the top-quark; all other quarks are coupled to the 
1st Higgs doublet (^i). This scenario is, therefore, 
realized by setting in (|23"|): 



Gij x (<5j-i + S j2 ) 
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G, 
. 



o j3 



(24) 



where G is again an unknown Yukawa 3x3 matrix 
in quark flavor space. 

Using' the Lagrangian pieces given above, we list in 
Table | all the couplings required for the calculation of 
T(Z -> didj) in a 2HDMII and in a T2HDM. 
In Table Q, sw> cw are the sine and cosine of the weak- 
mixing angle 9w, — Ta u , ,m t for i = 1,2,3, 
respectively, and 



Y — 



V2s w Mi 



— tan/3Vf 



CKM 



(25) 



with Vckm the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix and tan/3 = tp — v%lv\ (we will loosely refer to 
the ratio Vijv\ either as tan/3 or tp), where i>2(fi) is the 
VEV of &2{&i)- Also, the matrix £ introduced in Table 
| is composed out of the unitary matrix that diagonalizes 
the right-handed up-type quarks and the Yukawa matrix 
Ufj defined in (pif). It, therefore, arises from the specific 
structure of the Yukawa interactions in the T2HDM and 
can be parametrized as follows |l7j (neglecting the mass 
of the first generation up-quark): 



m r e 





m c e ct C 
mt (1 - \ecS 2 ) 



, (26) 



where e c t = m c /m t and £ is an unknown parameter (as- 
sumed here to be real) whose "natural" size is of 0(1). 

Notice that the specific structure of the T2HDM's 
Yukawa potential does not give rise to tree-level flavor 
changing couplings between a neutral Higgs and a pair of 
down-quarks (while allowing for tree-level neutral Higgs- 
top-charm couplings). Therefore, the decay Z — > didj is 
not affected at one-loop by flavor changing neutral Higgs- 
quark interactions. 



A. Z -> bs in 2HDMII 

The charged Higgs one-loop contribution to the decay 
Z — > bs in a 2HDMII was examined before in Q|. Here 
we wish to recapitulate the salient features of this decay. 

On the left side of Fig. || we plot BR(Z — > bs + bs) as a 
function of tan j3 for charged Higgs masses of 100, 400 and 
600 GeV. As can be seen, BR(Z — > bs + bs) is maximal 
for low tan/3 ~ 0(1) for which it is controlled by the top- 
quark Yukawa coupling which is oc 1/ tan (3. Thus, in this 
range the required flavor transition is mediated by t — * s 
and the BR(Z — ► bs + bs) is, therefore, essentially pro- 
portional to (?7Jt/ tan/3) 4 x (VQ b KM ) 2 (Vq KM ) 2 which is 
the square of the product of the tbH + and tsH + Yukawa 
couplings. 



At around tan/3 ~ 13 there is a "turning point" 
at which the BR(Z — ► bs + bs) starts to in- 
crease with tan/3. At this point the contributions 
from the top and charm-quark loop exchanges be- 
come comparable, since the charm-quark effect being 
oc m 2 m 2 tan 4 P{V£ s KM ) 2 {V£ b KM ) 2 (for tan/3 > 13 the 
charm-quark exchange is dominated by the Yukawa cou- 
plings b 2 ^ and bj^, see Table [j| equals that of the top- 
quark. As tan/3 is further increased (tan/3 > 13) both 
the top and charm-quark loop exchanges are dominated 
by the right-handed down-quark Yukawa couplings b 1 ^ in 
(||) (which is oc tan /3) and are, therefore, comparable and 
increasing with tan j3. 

Note that the curves in Fig. | for the 2HDMII sce- 
nario (the left side) pass through unrealistic values in 
the tan/3 — m H + plane. In particular, the decay b — » sj 
imposes strong constraints on the tan/3— m H + plane [p~8| : 
rriH+ ~ 400 GeV independent of tan/3. Thus, if tan/3 = 
1, then the largest allowed value for the BR(Z — > bs + bs) 
is ~ 10~ 10 if rrifj+ lies close to its lower bound from 
b — * 37. For even smaller values of tan (3, say tan f3 £ 0.5, 
the b — > S7 constraint requires m H + 
the BR(Z -> bs + bs) is again < 10 



> 500 GeV for which 
" 10 in the 2HDMII. 



B. 



6s in T2HDM 



The main difference between the T2HDM and the 
2HDMII charged Higgs sectors lies in the structure of the 
cdiH + Yukawa interactions (di — d, s, b for i = 1, 2, 3 
respectively). In particular, while in both models the 
top Yukawa coupling to down-quarks, i.e., the tndLH + 
coupling, is oc nitVcKMl tan/3, the charm-quark Yukawa 
coupling is completely altered by the presence of the 
matrix S in (|2^). Specifically, the Cfl6iiJ + coupling 
is oc m c ^*Vp^ M tan/3 and the crSlH + coupling is oc 
m cVc S KM tan/3. These couplings are, therefore, a fac- 
tor of tan 2 /3 x (Vq KM /V^ b KM ) and tan 2 /3, respectively, 
larger than in the 2HDMII if £ is of 0(1). 

As can be seen from Fig. |^ (the right side), in the range 
tan j3 £ 5 the BR(Z — > bs + bs) is practically identical in 
both the T2HDM and the 2HDMII; in this range it is 
dominated by the top-quark loop and, therefore, by the 
top-quark Yukawa couplings to the b and s-quarks which 
are essentially the same in these two versions of a 2HDM. 
On the other hand, for larger values of tan/3, in contrast 
to the case of a 2HDMII, in the T2HDM the charm-quark 
loop starts to dominate due to the enhancement in the 
CR,bLH + and crslH + Yukawa couplings (see discussion 
above). In fact, the cubLH + coupling is doubly enhanced 
in the T2HDM; first by the tan/3 factor and second by a 
factor of Vc KM /Vq KM , i.e., in this model this coupling- 
does not suffer from the CKM suppression factor Vckm ■ 

It should be emphasized that a large tan/3, e.g., 
tan/3 ^0(10), is the "working assumption" of the 
T2HDM. In particular, the T2HDM looses its attrac- 
tiveness in the small tan /3 regime, since in this range 
it fails to explain the large top-quark mass - this being 
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scalar (Soi—i) 
fermion (fi) 



~ l L(Zf) 
l R(Zf) 



°£(a=l) 



°R(a=l) 



a=ld=l 

9z 



2HDMII 



H+ 
Ui, i = 1,2, 3 

a-L(Zu) 
a-R(Zu) 



% 3 " tan^ 3 " 



y^7 x 



T2HDM 



m±Vij 



H+ 

Ui, i = 1,2, 3 

1fi(Z«) 

1 + 



ifc V CKM | T _i_ 1 

tan 7 /3 



TABLE I: The couplings required for the calculation of T(Z -> did./) in a 2HDMII and a T2HDM. The couplings follow the 
notation used in the Feynman rules of (^), (Q) and (^). Also, a fl ( Ztl ), a L ^ Zu ^ are given in (|^)-(|^). 

2HDMII T2HDM 
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10 
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111,,' LOO GeV 
m H *=400 GeV 
III,, (>(K) GeV 




M i 



O IO 20 30 40 50 O IO 20 30 40 50 

f ; 1 1 ■ | > 

FIG. 2: -> 6s + Is) as a function of tan/3, for m H+ = 100, 400 and 600 GeV, in a 2HDMII (left side) and a T2HDM 

with £ = 1 (right side). 



the main motivation behind this version of a 2HDM. At 
the same time, taking into account low-energy experi- 
mental data from K — K mixing, ex and b — > sj, the 
tan (3 — m H + plane is also constrained in the T2HDM 
jl7], [l9) , especially in the large tan (3 region in which this 
model differs from the usual 2HDMII. For example, for 
£ = 1 and taking the SM best fit value for the Wolfen- 
stein parameters p and 77, then the e^- constraint implies 



m H + > 500 GeV for tan/3 - 20 and m H + > 4 TeV for 
tan f3 r~j 50 |L9|_. Imposing these bounds we find that 
BR(Z -> bs + bs) ~ O(10 -8 ) is the best case value in 
this model assuming a large tan /3. 
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IV. SUPERSYMMETRY WITH SQUARK 
MIXINGS 



soft breaking Lagrangian terms in the squark sector: 



In SUSY flavor changing phenomena can emanate from 
mixings of squarks of different generations through the 



J 



squark 
soft 



I 



(27) 



where Q is the SU(2) scalar doublet and U, D are the up- 
squark, down-squark SU(2) singlets, respectively. Also, 
i,j are generation indices. 

The mass matrices in the up-squark and down-squark 
sectors may then be generically expressed as: 



M- - = 

U,D 




«n)L 2 



(28) 



where (m^) LL , 



(m| q)rr and (m± q)lr are 3 x 3 ma- 
trices in squark flavor space. In the super CKM basis 
the squark fields are rotated "parallel" to their fermionic 
super-partners. In this basis, and assuming that there is 
a typical common mass scale for the squarks, mo, which 
is sufficiently heavier than the electroweak mass scale 
(mg/Af| 3> 1), these matrices are related to the soft 
breaking bilinear and trilinear terms in (E7p via p(| : 



= V?M%Vp 



( m fr)LL 



- V%M? T T V R 



r 2n _ "sin/3 ^ A *„uj 
{™>f r )LR = —VI u A„,V„ 



V2 



u y R 



(ml) LL = V L D M 2 Q V^ , 



K) 



LR 



(29) 



r 



where V L D R and V^ R are the rotation matrices that di- 
agonalize the down and up-quark fermion mass matrices, 
respectively (the CMK matrix is Vckm = VE^l^)- 

Assuming that flavor changing squark mixings are sig- 
nificant only in transitions between the third and second 
generation squarks, we choose the following textures for 
the 3 x 3 matrices (m| ^)ll and (m| q)rr: 




1 

j[/,Z>(32) 
J LL,RR 





M,D(23) 
LL,RR 

1 



l ) 



(30) 



where 5 



(7(23) M(32) M(23) 



,D(23) 
U RR ' 



LL 

,D(32) 



'LL 



RR 



M(32) 
RR 



and S 



D(23) 
LL 



,D(32) 
J LL 



5 RR represent flavor mixings in the t — c and 



b — s sectors, respectively. These flavor violating quan- 
tities emanate from non-diagonal entries in the bilinear 
soft breaking terms Mq, M^j and Mf, in (|27|). Simi- 

are associated with 



larly, 5 U L { R \ 5 U L ( R 2) and <5 i? 
non-diagonal (flavor changing) entries in the trilinear soft 



D(23) ^73(32) 



LR. 



breaking terms A u and Ad defined in ( |27| ) . As 
adopt the following simplified Ansatz: 



VZKVr 




and 




111 we 



(31) 



(32) 



such that A is a common trilinear soft breaking param- 
eter for both up and down-squarks and the parameters 
XutVu and Xd,yd represent flavor mixing effects in the 
t — c and b — s sectors, respectively. It then follows that 
d LR , o LR and 6 T i, , 6 T i, are related to x,,,, y„, and 
x d , yd via: 



(33) 



,D(23) 
} LR 



COS/3 vA £,(32) 

x d r— X 9 ' " T 

V2 "Iq 



cos (3 vA 



8 

(34) 



Within this mixing scenario, in which flavor changing 
effects in the squark sector are present only in the second 
and third generations, the 6x6 mass matrices in the up 
and down-squark sectors reduce to 4 x 4 matrices. 

For the t — c sector, in the basis = (cl, c r , t L , t R ), 
we then have: 



M 



( 1 








1 


,1/(32) 
LL 

Wr 


,17(32) 

LR 
,17(32) 

RR 



.17(23) 


,17(23) 


LL 


°LR 


,(7(32) 


,17(23) 


°LR 


°RR 


1 


-X t /m 


-X t /ml 


1 



/?? 



i 



(35) 



and similarly for the b — s sector, in the basis $^ 
(sl, Sr, b L , b R ), we have: 



( 1 



Mi = 



5 





73(32) 
LL 
D(23) 
LR 




1 

A r>(32) 

°LR 
,77(32) 

RR 



,£>(23) 
LL 

,D(32) 
LR. 

1 



,15(23) \ 

,D(23) 
RR. 

-X b /ml 
1 / 



(36) 



The factors and Xf, are responsible for mixings be- 
tween left and right handed squarks of the same genera- 
tion and are given by: 



X t = V -^A 



X h = - 



(37) 

A + nib ta,n [3 fi , (38) 



V2 tan/3 ' 

v cos 



where n, is the usual Higgs mass parameter in the SUSY 
superpotential. 

After diagonalization of M%t an d ^sb one obtains the 
new mass-eigenstates which are now t — c and b — s ad- 
mixtures, respectively. The diagonalizing matrices Rjj 
and Re, arc defined through: 



R-D,ik®D,k 



(39) 



where 




and ul,r, <1l,r (u = c,t and d = s,b) are the SU(2) weak 
states, while ui,2, di,2 are the corresponding physical 
states (mass-eigenstates) . 

Let us now consider separately the cases in which the 
one-loop flavor changing decay Z — > bs is driven either 
by the t — c or by the b — s mixing phenomena. 



A. b — S mixing 

Here the flavor changing decay Z — > bs is generated 
by one-loop exchanges of the b— s admixture states, &d, 
and gluinos, g. We thus have S a = &D,a with a = 1 — 4, 
and / = g in the diagrams of Fig. [I]. Note that diagram 
(2) in Fig. [I], which requires the Vff coupling, does not 
contribute since the Z-boson does not couple to gluinos 
at tree- level. 

The one-loop b — s mixing effect on the decay Z — > bs 
was considered more than a decade ago in @], where it 



was assumed that flavor violation in the down-squark sec- 
tor is controlled by radiative corrections to the down- 
squark mass matrix induced by off-diagonal CKM ele- 
ments. Instead, as described above, we assume here that 
the flavor violation is rooted with non-diagonal entries 
in the soft SUSY breaking sector. The approach taken 
here is, therefore, fundamentally different from the one 
suggested in 

The following interaction vertices are required for the 
calculation of Y(Z — ► bs) in the b — s mixing scenario: 

V^S a S p -» Z^I^dj , (41) 
These are derived from [p2[: 



[1] Note also that B used unrealistic top-quark and squark masses. 
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A u H B M 

2 Svy 6 Cw 

1 e - «"t - 
-i- — BpdRj d'i d* Rl , 

6 Cyi/ a 



dt, m d L , 



(42) 



C{dgd) = g s V2T a g (-d* L<e L + d* Re Rj dj + h.c. , (43) 

where g s is the SU(3) coupling constant and T a are the 
SU(3) group generators. 

Using the above interaction Lagrangian terms, the cou- 
plings required for the calculation of T(Z — > 6s) in the 
form defined in (0) and fll) are obtained by rotating the 
weak states, to the physical states, <&£>, according to 
(|39|). These couplings are given in Table |§ 



SUSY with b — s mixing 



scalar (S a ) 
fermion (fi 



1 L(Zf) 
1 R(Zf) 



b 1] 



b ZJ 

°R{o.) 



9? 



-D,a, 



-V2g s T a {R* D>3a 5 3j + J£, la 4y) 
-V2g s T a {Rh Aa h 3 + R* D , 



TABLE II: The couplings required for the calculation of 
T{Z -> bs) in the MSSM with b-s mixing. These couplings 
follow the notation used in the Feynman rules defined by (m) , 
(§) and (§). 



The relevant low-energy SUSY parameter space for the 
6 — s mixing scenario is characterized as follows: ^ 

I. Some of the flavor changing parameters in the 6 — s 
sector are severely constrained by the 6 — > sj decay 
@|2|. In particular, <5^ 23) , 6%£ 2) < O(10~ 2 ) is 
required in order to keep BR(b — > S7) within its 
experimental measured value. I 3 ' On the other hand, 



C 3) -d ST' ofO(l) are not 
ruled out by 6 — > S7 nor by any other low energy 
process that we know of. 

We will, thus, use the LL and RR delta's as the 
only source for b — s mixing. Moreover, since there 
is no a-priory theoretical reason for the four LL and 
RR delta's to be significantly different, we will set 
all of them to a common value denoted by 5 D . That 



is, we fix 6 



D(23) 



? r>(32) 



xD(23) _ 

and vary 6 in the range < 5 < 1. 



•D(32) 

'rr 



II. The SUSY parameter space needed to evaluate 
T(Z — > 6s) in this scenario is: mo, /Li, ^4, tan/3, 
and J D . The low-energy values of these six param- 
eters fully determine the 6 — s scalar spectrum (i.e., 
masses and mixing matrices) and the gluino mass 
(m§), from which all the couplings in Table || are 
calculated. These six parameters will be varied sub- 
ject to the requirement that squark masses as well 
as the gluino mass are heavier than 150 GeV. 

In Figs. |, I and @ we plot BR{Z -> 6s + 6s) as a 
function /Lt, S and tan/3, respectively, for three values 
of the common squark mass: too = 1000, 1600 or 2200 
GeV. The rest of the parameters are varied subject to 
the above criteria |2J] . In order to better understand the 
dependence of BR(Z — > 6s + 6s) on the physical squark 
mass spectrum, we accompany Figs. || and |^ by Figs. |] 
and ||, respectively, in which we depict the masses of the 
four physical squarks ms 1 2 and that correspond 

to the same choices of the SUSY parameter space as in 
Figs. I and |. 

Let us summarize the results shown in Figs. |3U 



• The branching ratio of the decay Z — * 6s is en- 
hanced dramatically with the increase of the mass 
splittings between the four physical squarks. This 
is due to a GIM-like cancellation which is opera- 
tional in the limit of degenerate squark masses as a 
result of the unitarity of the rotation matrix -RdJ 4 ' 
Thus for example, a typical mass spectrum that can 
drive the branching ratio to the 10~ 6 level is when 
the lightest down squark, b\, has a mass below 250 
GeV, while the rest of the squarks have masses at 
the 1-3 TeV range. 

• As expected, BR(Z — > 6s) drops sharply as S D is 
decreased. Clearly, this is traced to the fact that 
the mixing among the bottom and strange type 
scalars diminishes in the limit S D — ► 0, see Figs. || 
and ||. 



[2] The term "low-energy" refers to the electroweak (or collider ener- 
gies) scale and is used in order to distinguish it from the scale in 
which the soft breaking couplings are generated (e.g., the GUT 
scale). 

[3] The bounds on the different delta's reported in pCl, E3| were ob- 
tained using the mass insertion approximation, wmle we perform 
an exact diagonalization of the squark mass matrices. There- 



fore, in the cases where 0(1) delta's are allowed (e.g., for S L ^ ) 
these bounds may only serve as indicative for their expected size, 
since the mass insertion approximation necessarily assumes small 
delta's. 

[4] The unitarity of Ro also ensures that the infinities that appear 
in the individual diagrams of Fig. hi cancel. 
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SO , A=lSOO GeV 




FIG. 3: BR(Z — > bs + bs) as a function of the Higgs mass parameter /i, for combinations of mo = 1000, 1600, 2200 GeV with 
m- g = 200, 600 GeV, for tp = 50, A = 1500 GeV and for 8 D = 0.5 (left plots) and 8 D = 0.9 (right plots). 



S"=o.s» 
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FIG. 4: Physical masses of the second and third generation down-type squarks as a function fi. The rest of the parameters are 
as in Fig. H. 



• For a sufficiently large tan/3, BR(Z — > 6s) is al- 
most insensitive to the value of the common trilin- 
ear soft breaking parameter A as long as (j, is large 
enough to drive the desired mass splittings between 
the squarks. This behavior is due to the dominance 



of the /i term in the quantity Xb defined in (|3^) for 
large tp (recall that is responsible for the mass 
splitting between the two bottom- type scalars). On 
the other hand, for tg ~ 0(1), the term tx A in Xb 
(see the r.h.s. of (|38|)) becomes important when 




FIG. 5: BR(Z — > bs + bs) as a function of the flavor mixing parameter 5 D , for combinations of mo = 1000, 1600, 2200 GeV 
with /i = -1000, - 2000, - 3000 GeV, for tp = 50, A = 1500 GeV and for m g = 200 GeV (left plots) and m- a = 600 GeV 
(right plots). 
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FIG. 6: Physical masses of the second and third generation down-type squarks, as a function of S D , for mo = 1000, 1600, 2200 
GeV with n — —1000 GeV (left plots) and /i = —3000 GeV (right plots). The rest of the parameters are as in Fig. ^. 



A ~ fi. This feature can be seen in Fig. 0. 

• For the reason outlined above, BR(Z — » bs) is sym- 
metric about /i = for large tan/?, in which case 
the term oc jj, in X}, dominates and the effect of the 



A term is negligible. 

• For p?jA 2 >> 1, BR(Z — > 6s) is increased with 
tan/3. Again, this is related to the dominance of 
the fj, term in for large tan/3. 



A=1000 t ; t-V A=3000 GeV 

S D <>_•> , m==200 GeV 
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FIG. 7: -> &s+6s) as a function of tan /3, for combinations of m = 1000, 1600, 2200 GeV with /i = -1000, -2000, -3000 

GeV, for m s = 200 GeV, 8 D = 0.9 GeV and for A = 1000 GeV (left plots) and A = 3000 GeV (right plots). 



• BR{Z — > bs) drops with nig. 

To conclude this section, we have shown that BR(Z — > 
bs) ~ 0(1O~ 6 ) can be achieved in the b — s mixing 
scenario provided that the gluino mass and one of the 
third generation down-type scalar masses lie close to the 
electroweak mass scale, while the rest of the down-type 
squark masses are at the TeV range, i.e., a large mass 
splitting between the lightest and rest of the down-type 
squarks is needed. Such a mass hierarchy in the squark 
sector requires a typical "heavy" SUSY mass scale with 
soft breaking parameters at the level of a few TeV. This 
scenario is somewhat motivated by the non-observability 
of SUSY particles in past and present high energy collid- 
ers. 



the t — c admixture states, <&[/, and charginos, \- More 
specifically, we have S a — <&u,u with a = 1 — 4, and 
fi = Xi with i = 1,2 for the two charginos (we find 
it convenient to calculate the exchanges of the charged 
conjugate chargino states x c )- 

Thus, in the t — c mixing scenario the following inter- 
action vertices are required: 



B. t-c mixing Vnfifj ~> Z nXiX C j 



V„S a Sp -> Z^®*^ , (44) 



In the stop-scharm mixing scenario the flavor changing S a fidj — > ^qqXMj 

decay Z — ► bs proceeds through one-loop exchanges of These vertices are taken from p2f: 



c(v^x-x c j) 

C{V^uu) 
C(u X c d) 



X C jl^{a[ {VxC) L 



2sw c w 



= — % 



2 sw 



UL.ed 3 (/ 



6 cw 



«-» 2 e M 

o cw 



Htij) T _l t R ( a 3)fi\ x C VcKM 



+ UR. 



r 



(45) 
(46) 
(47) 



where 



' l L{Z X ") 



Z~Z-* + cos29 w S, 



(48) 
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SUSY with t — c mixing 



scalar (S a ) 
fermion (fi 



$r/,a, a = 1,2,3,4 
Xi, i = l,2 



1 J 

2 L(Zf) 
1 R(Zf) 



ti 3 



b 13 



ol8 

9z 



»3 



2s w c w a HZ x ") 
2s w c w u R(Z x c ) 

f L ' ^ {Ru,lc,VcKM + Ru,3a ^OKm) + { m <:Ru,2ci^CKM + m tRu,ia^( 



33 ) 
CRM) 



TABLE III: The couplings required for the calculation of Y{Z —* bs) in the t — c mixing scenario. The couplings follow the 
notation in (^), (]?]) and Also, a^( ZxC -. and a^ z ^ are given in ( fl8| ) and (|49|), f^ m ^\ f^ m ^ and f^ < ~ al3 ' > are defined in 
(j5(i|). is the rotation matrix defined in ([39|), 



<»_•> . 5 l 



< > 



o , S l 



l <> - 
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FIG. 8: B_R(Z -> &s + 6s) as a function of tan /3, for combinations of m = 1000, 1600, 2200 GeV with m 2 = 200, 600, 1000 
GeV, for A = 1000 GeV, fx = -2000 GeV and for 5^ = 0.9, 8% = (left plots) and = 0, 5% = 0.9 (right plots). 



4W) = (^i +cos20 w %) , 



f L(lij) 
J L 

f R(Uj) 
J L 



m dj 



y/2s w M w cos (3 



J 2i i 



SW 



(49) 



JR. 

MHj) 

JR 



y/2s w M w sm[3 



(50) 



and are the chargino mixing matrices given in [ p2[ . 
Also, A 3 and 5 are the SU(2) and U(l) gauge fields, 
respectively, and Ul,r are the SU(2) weak states of the 





squarks 




m,,=1000 GeV 






■ ■ l a ' 




1 1 1 i , 




m r 


m„ I <><><> GeV 








m„=2200 GeV 











tliai<iiiios 
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FIG. 9: Physical masses of the second and third generation up-type squarks (left plots) and of the charginos (right plots), as a 
function of tan (3. Squark masses are shown for mo = 1000, 1600 and 2200 GeV and for 5m — 0.9, 5^=0 (equivalent to the 
case 8 M = 0, <5^ = 0.9, see text) and chargino masses are given for mi = 200, 600 and 1000 GeV. The rest of the parameters 
are as in Fig. 



up- type scalars. 

Here also, the couplings needed for the calculation of 
T(Z — > bs) in the form defined in (||), (0) and (||) are 
obtained from the Lagrangian terms in (|45|)-(|50|) by ro- 
tating the weak states, to the physical states, 



accor ding to 
Table Ell. 



I). These couplings are summarized in 



The contribution of the t — c mixed states to the one- 
loop diagrams in Fig. [I] are characterized as follows: 

I. The quantities that mediate the flavor changing tran- 



sition b — > s in the t — c mixed scenario are: 5 



,17(32) 
J LL 



,17(23) 
RR ■ 



,(7(32) 
RR 



,1/(23) 
J °LR 



and 



U(23) 
LL ' 



Recall that 



the LL and RR delta's originate from the bilinear 
soft terms in (|27]), while the LR delta's are associ- 
ated with the trilinear soft breaking SUSY terms. 
Thus, we will separate these two types of flavor 
violating sources in our numerical analysis. In par- 



ticular, we define 5V* — 5 



U(23) 



J LL 



X U(23) _ 



cU(32) 
LR 



,U(23) 
RR 



,17(32) 
°LL 

LR — u LR and we vary ei- 
ther 5m or 5 a in the range [0,1]. Note that an 
0(1) value for either 5jL or 5% is consistent with 
all experimental data [EO, p3|. 



II. 



The required SUSY parameter 
mo, /i, A, tan/3, m2 and 5^, 



space is: 
5%> where 



ni2 is the SU(2) gaugino mass parameter. The 
low-energy values of these six parameters fix the 
t — c scalar spectrum (i.e., masses and mixing 
matrices) and the chargino masses and mixing 



matrices from which all couplings in Table [II are 
derived. 

As in the b—s mixing case, these parameters will be 
varied subject to the requirement that the squark 
masses are heavier than 150 GeV and, in addition, 
that the charginos are heavier than 100 GeV p4fl. 



Taking maximal flavor violation in the t — c mixing 
scenario, i.e., 5^ ~ 0(1) or 5% ~ 0(1): an d varying 
the rest of the SUSY parameters involved subject to the 
above criteria, we find that BR{Z — * 6s + bs) can reach 
fewx 10 -8 at best. Here also, the BR(Z — > bs + bs) is sig- 
nificantly enhanced when large mass splittings between 
the four up-type squarks, 2 and mj are present. 
Such a hierarchy in the up-type squark mass spectrum 
makes the GIM-like cancellation mentioned earlier less 
effective. 

Indeed a two orders of magnitude difference between 
the t — c and b—s mixing cases is expected due to an 
(a s /a) 2 enhancement factor in the b — s scenario (com- 
pared to the t — c mixing case) which arises from the 
gluino QCD coupling. 

In Figs. H we plot BR(Z — > bs + bs) as a function 
of tan/3, for combinations of mo — 1000, 1600 and 2200 
GeV with m 2 = 200, 600 and 1000 GeV and for either 
{<% = 0.9, 5% = 0} or {8% = 0, 5% = 0.9}. [5] For illus- 



[5] Note that the physical up-squark masses have the same depen- 
dence on 8 M or S'i when one of the two delta's is set to zero. 
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tration we set A = 1000 GeV and fi = -2000 GeV. In 
Figs. ^ we depict the masses of the four physical up- 
type squarks 2 and and the masses of the two 
chargino states as a function of tan /3, for the same SUSY 
parameter choices as in Figs. || Q. 

V. SUSY WITH R P VIOLATING 
INTERACTIONS 

If Rp is violated in the SUSY superpotential, then fla- 
vor changing transitions can emerge from interactions of 
squarks or sleptons with fcrmions. In particular, there 
are two types of RPV terms that are allowed in the su- 
perpotential if the discrete Rp symmetry is not imposed. 
These are the RPV trilinear Yukawa-like (RPVT) oper- 
ators and bilinear (RPVB) operators. 

In the usual convention, the RPVT are proportional 
to the dimensionless couplings A, A' and A", see e.g., 
p5[ . Here we will assume that A" <C A' and investigate 
the one-loop effects of the A' type operator on our flavor 
changing decay Z — + djdj:^ 

Wrpvt D e ab \' ijk L?Q b j D c k , (51) 

where Q and L are SU(2) doublet quark and lepton su- 
permultiplets, respectively, and D is the SU(2) singlet 
down-type quark supermultiplet. Also, k = 1,2 or 3 
are generation indices and a, b are SU(2) indices. 
The RPVB operator is: 

Wrpvb = -e abl i t L?H b u , (52) 

where H u is the up-type Higgs supermultiplet and i — 
1, 2 or 3 labels the lepton generation. 

In addition, if one does not impose Rp, then the usual 
set of Rp conserving (RPC) soft SUSY breaking terms is 
extended by new trilinear and bilinear soft terms which 
correspond to the RPV terms of the superpotential, i.e., 
to the ones in ([[l]) and ([52"|). For our purpose, only the 
followingsoft SUSY breaking bilinear term is relevant 
© 13,11: 

Vrpvb = ^abbiL^H^ , (53) 

where L and H u are the scalar components of L and H u , 
respectively. 

The RPVT operator (oc A') in ( pl| ) gives rise to the 
following scalar-fermion-fermion RPV interactions: 



c = Kjk {*i44 + 444 + (4)* (4) c 4 

-§14< - 444 - (4) * (el) C <} + h.c. . (54) 

where d(u) is a down(up)-quark, e(y) is a charged- 
lepton(neutrino) and scalars are denoted with a tilde. 

The RPVB operator (oc /i,) in ([52]) gives rise to mix- 
ings among charged leptons and charginos as well as be- 
tween neutrinos and neutralinos. However, low energy 
flavor changing processes [p9| , flavor chan ging leptonic 
Z-decays |}0| and neutrino masses |2(| |3(], |3l|, |3^] sug- 
gest that the pn are expected to be vanishingly small. 
We will, therefore, neglect its contribution to the decay 
Z — > djdj\ 7 ^ On the other hand, the soft breaking RPVB 
term (cx bi) in (^3|) gives rise to mixings between sleptons 
and Higgs-bosons which may be exchanged in the loops 
of the diagrams shown in Fig|j]. 

Let us now categorize the different types of RPV inter- 
actions that contribute at one-loop to the flavor changing 
decay of Z — > didj. Since the decay Z — > didj conserves 
Rp, there should be two insertions of RPV vertices in 
the one-loop diagrams of Fig.[j]. We can thus divide the 
various types of RPV one-loop exchanges into two cate- 
gories, type A and type B, according to the pair of RPV 
couplings involved: 

Type A: The RPV contributions that are proportional 
to the product A'A', i.e., T(Z — > didj) oc A'A', 
where A' is defined in (|5l|). 

Type B: The RPV contributions that are proportional 
to the product 6A', i.e., T(Z — > didj) oc b\' , where 
b is the soft breaking RPV bilinear coupling defined 




A. Type A RPV effect 

The type A RPV contribution to Z — * djdj emanates 
from the first five RPV Yukawa-like interaction vertices 
in (|54|). In this case we assume that 6,; — > such that 
mixing effects between sleptons and the Higgs fields are 
absent. 

We can further sub-divide the type A contributions 
into 6 types according to the type of scalar (5) and type 
of fermion (/) that are being exchanged in the loops: 

type Al : S a = e L , a , fi = u t 
type A2 : S a = d L . a , f. L = Wj 



Thus, for example m fii (5^ f = 0.9, 5^ = 0) = m a .(<5^- =0,S% = 
0.9), i = 1,2,3 and 4. * 
[6] Note that at the one-loop level the A type couplings do not con- 
tribute to the decay Z — > djdj. 



[7] The one-loop exchanges of possible lepton-chargino and neutrino- 
neutralino admixture states in Z — > djdj will be controlled by 
the square of the RPV couplings product \Xi X A'. 
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type A3 : 


S a 


— d,R ia 


fi = ^ 


type A4 : 


S a 


= VL.a , 


fi = di 


type A5 : 


S a 


= v'l, a , 


fi = di 


type A6 : 


S a 




fi — &i 



(55) 



where a = 1,2, 3 and i = 1, 2, 3. 

For each of the type A RPV exchanges above, the 
generic couplings defined in (||) , (]?]) and (|J) are summa- 



rized in Table 



In particular, for a given /, the Zff 
couplings of (|2|) are given by (§)-@- The Sdf couplings 
are taken from the Yukawa like interactions in ( |54| ) , while 
the ZSS couplings are extracted from CiV^uu) in (46), 
from £(V^G?d) in p^ ) and from the VLL interaction La- 
grangian: 



c(yJLl) 



— T 3 A 3 - 

~ A a "II 

Sw cw 



-B„ 



(56) 



where A 3 and B are the SU(2) and U(l) gauge fields, 



respectively, L 



Given the couplings in Table |lVj and the structure of 
the form factors in 



fllOD-fllq) it is evident that there are 
only two types of A' A' product combinations which enter 
the type A RPV contribution to the decay Z — ► didj: 

1. The product A^ nJ A^ nJ . Types Al, A2, A5 and A6 
are proportional to this couplings product. 



2. The product X'^ lIn X' mJn . Types A3 and A4 are pro- 
portional to this couplings product. 



Furthermore, since none of the scalars have both a 
left and a right handed RPVT coupling to fermions in 
the type A scenario, i.e., in the notation of (||) either 
^l(q) — or ^R(a) = ( see Table IV), the form fac- 
tors Bjf k and Bj^ k in the amplitude (|9|) (which requires 
a none-zero value for both the left and the right handed 
scalar-fermion-down quark couplings) vanish. Also, since 
= for the RPV contributions of types Al, A2, A5 
and A6, they contribute only to the right-handed vector- 

the RPV contributions 

0, therefore, contribut- 



like form factor A^ k . Similarly, the RPV contributions 
of types A3 and A4 have ti^rs 
ing only to A^ k . 

It should be noted that for any one of the type A RPV 
exchanges, if the scalars of different flavors that are be- 
ing exchanged in the loops are degenerate and upon ne- 
glecting all fermion masses except for the top-quark, then 
there remain only three distinct types of contributions of 
the A' products in the type A RPV scenario. That is, un- 
der this assumption BR(Z — ► didj) can have only three 
different values which we denote by BRl IJ , BR2 IJ and 
BR3 IJ as follows: 



BR1 IJ = BR(Z -> didj) when (A' y/ x X'^jf ^ ; j ^ 3, % = 1,2,3 (57) 
BR2 IJ = BR{Z djdj) when (A- 3J x \' aJ f ± ; i = 1, 2, 3 (58) 
BR3 IJ = BR(Z -> didj) when (A-^ x X-^) 2 ^ ; i,j = 1,2,3 , (59) 



such that BR{Z -> bs) = BR1 32 , BR2 32 or BR3 32 de- 
pending on which of the three A' A' product combinations 
is non-zero. 

In Table ^ we give a sample of our numerical results 
for the three BR's in (|57[)-(^9[) scaled by the square of 
the appropriate A'A' product. The results presented in 
Table [V] correspond to the case of a single non-zero A'A' 
product (one index combination) contributing to each of 
the BR's BR1 32 , BR2 32 and BR3 32 . In addition, the 
masses of the squark and slepton being exchange in the 
loop (for a given index combination of the corresponding 
A'A' product) are set to either nig = 500 GcV with m» = 
200 GeV or m§ = 1000 GeV with m E = 500 GeV. 

The existing limits on the A' coupling products in fl57|)- 
(|S9|) seem to indicate that the typical allowed values of 
A' x A' for any of the index combinations in (|57|)-(|59]) is 
at the level of ~ few x 10~ 2 It should be noted, 

however, that the limits reported in |33j assume 100 



GeV scalar masses. These limits scale with the scalar 
masses (typically as [m^/100 GeV] 2 , where is the 
scalar mass) and are, therefore, relaxed as the scalars 
become heavier. ^ 

Using A' x A' ~ O(10~ 2 ) in conjunction with the results 
presented in Table 0, we see that the expected branching 
ratio for Z — > bs in the type A RPV scenario investigated 
in this section lies in the range BR(Z — > djdj) ~ 10~ n — 

10- 10 . 

This type A RPV one-loop effect in BR(Z ->■ didj) 
was also investigated in || . Although || evaluated some 
distinct limiting cases of the type A RPV contributions, 



[8] Note that b — > sy, which is proportional to A' X A' products with 
the same index combinations as in Z — > bs, allows some of the 
above A' X A' coupling products to be at the 10 — 1 level [p4|. 
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type Al 


type A2 


type A3 


scalar (S a ) 


&L,a, a =1,2 


dz-.a, a = 1, 2, 3 


dfl.a, a = 1, 2, 3 


fermion (fi) 


Mi, i = 1, 2, 3 


i>i, i = 1, 2, 3 


i/f, i = 1,2,3 


ij 

a L(Zf) 






—OiR{Zv) 


U R(Zf) 






— 0-L{Zv) 


b ij 








v 

/ »ija 




y* 







9z 


-efr-'WSci, 




3 c^i/ ^ 


2s\yC\y 







type A4 


type A5 


type A6 


scalar (S a ) 


z^a, a = 1,2 


vl t a, a = 1,2 


a = 1, 2, 3 


fermion (/;) 


di, i = 1,2,3 


di, i = 1, 2, 3 


e», i = 1, 2, 3 


a ij 

L(Zf) 




0-L{Zd) 




R(Zf) 




O-R(Zd) 




°L(c) 


V 








°H(a) 





V* 


,/* 


a/3 


e 


e ^ 






2s w c w ° a 


2s w c w i 1 3 S w) <W 



TABLE IV: The couplings required for the calculation of T(Z — > d/dj) in the type A RPV scenario. The couplings follow the 
notation in 







BR2 32 


BR3 32 




( A i33 xA i32) 




mg = 500 GeV, m £ - = 200 GeV 
mg = 1000 GeV, m- t = 500 GeV 


4.2 x 10" 7 
3.9 x 10" 7 


2.4 x 10" 6 
6.4 x 10 -8 


3.4 x 10~ 6 
3.0 x 10" 6 



TABLE V: Results for the three types of branching ratios BRl 32 , BR2 32 , BR3 32 as defined in each scaled by its 

appropriate A' A' coupling product. Results are given for two sets of squark and slepton masses as indicated. 



our results agree with the highlights of their analysis, i.e., perimentally allowed range of values for its relevant RPV 
that the typical BR(Z — > djdj) is expected to be at the parameter space, 
level of 1CT 11 - 1CT 10 if A' x A' ~ 0(\Q- 2 ). 

Thus, the type A RPV scenario is expected to yield a 
BR smaller even from the SM one. We, therefore, pro- 
ceed below to the second RPV type B scenario which 
seems to give a much larger BR(Z — > bs) within the ex- 
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B. Type B RPV effect 

The type B RPV effect arises when a Higgs particle 
that is being exchanged in the loops mixes with a slepton 
through the RPVB operator in (S) and then couples to 
the external down quark via a A^ype coupling. 

For simplicity we will assume that bi 7^ only for 
i = 3 in (|53|), thus, considering only the mixing between 
the third generation sleptons (L 3 ) and the Higgs scalar 
fields (Hd and H u )\^ It should be noted that b 3 7^ 
leads in general to a non-vanishing tau-sncutrino VEV, 
v 3 . However, since lepton number is not a conserved 
quantum number in this scenario, the H d and L 3 super- 
fields lose their identity and can be rotated to a particu- 
in which either /i 3 or v 3 are set to zero 
2l |35| . In what follows, we find it convenient to 
choose the "no VEV" basis, v 3 = 0, which simplifies our 
analysis below. 

Let us define the SU(2) components of the up-type 



lar basis (H 1 ,,, L' 3 



Higgs, down-type Higgs and L 3 scalar doublet fields (set- 
ting v 3 = 0): 



U 
,~,o ,-,0 



K 

i£ + v u +i<j>Z)/V2 



(60) 



where vH_, v 



and eg are the SU(2) CP-even, CP-odd 
r-sneutrino and left handed stau fields, respectively. 

When b 3 7^ the 3rd generation slepton SU(2) fields 
in j60| ) mix with the Higgs fields. In particular, in the 
basis $p = (h* , ,e$), the squared mass matrix in the 
charged scalar sector becomes: I 10 ! 



M 2 C = 



f c l[( m 

V 63 



a) 

• 



+ m 2 w ] s p c f3 [(m A ) 2 + 



nr 

+ m 



W-'J 



[( 



b 3 tp 



b 3 t 



■3 1/3 



(61) 



Mr 



cos 2/3, 



where m A and m sv are the pseudo-scalar Higgs mass and 
the tau-sneutrino mass, respectively, in the RPC limit 
63 - 0. 



Similarly, in the basis $° E = v+), the CP-even 

neutral scalar squared mass matrix becomes (at tree- 
level): I 1 



M 2 



(m A )-s 

[KO 2 



2 2 

I] s /3 c /3 



M/3 



[W 

K) 2 4 + 



ht/3 
-b 3 
(m° su ) 2 



(62) 



Finally, in the CP-odd neutral scalar sector and in the 
basis $q = ((p d , <j>®, ) one obtains: 



M 2 = 



{m A ) 2 c 2 p {m A ) 2 c p sp b 3 t{3 
{m A ) 2 cp S {3 (m A ) 2 sj 63 
b 3 



(63) 



b 3 tf3 



(ml 



The new charged scalar and CP-even and CP-odd neutral 
scalar mass-eigenstates (i.e., the physical states) are then 
derived by diagonalizing Mq, M| and Mq, respectively. 
Let us denote the physical states by: 



The consequences of b\ 7^ and/or 62 7^ is to introduce addi- 



$f7 = 




$ K = 









(i) 







(64) 



tional mixings among sleptons of different generations and mix- 
ings between the selectron and/or smuon scalar doublets with 
the Higgs fields. These extra mixing effects are not crucial for 
the main outcome of this section. 
[fO] We neglect the mixing between the right-handed SU(2) stau sin- 
glet and the charged Higgs fields which is proportional to the tau 
mass. 

[If] The one loop corrections to the 2x2 Higgs block in Mj, can 
cause a significant deviation to the tree-level mass of the light 
Higgs. This effect will be discussed below. 
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such that, for a small RPVB in the scalar potential, the 
new physical states in (|&4| ) are the states dominated by 
what would be the corresponding physical states in the 
RPC limit, 63 — 0, for which the Higgs sector decouples 
from the slepton sector in j6lj), (p3) and (|63|). In partic- 
ular, if 63 — * 0, then H, h, A and H + become the usual 
RPC MSSM's CP-even heavy Higgs, CP-even light Higgs, 
CP-odd pseudo-scalar Higgs and charged Higgs states, 
respectively. Similarly, in this limit and v T _ become 
the two mass-degenerate CP-even and CP-odd sneutrino 



states with a common mass 



= ™0 



while 



For the type Bl RPV contribution in ( pq ) the form fac- 
tors defined in (S) are calculated following the prescrip- 
tion described in section II. The generic couplings defined 
in (0), (Q) and (||) are summarized in Table VI for the 
type Bl RPV exchanges. In particular, the Sdf couplings 
(for S = <&c and f — u) are a combination of the Yukawa- 
like trilinear RPV interactions in (54) (those with the 
third generation slepton indices) and the charged Higgs 
Yukawa couplings which are the same as in the 2HDM of 
type II (given in section III). 

The ZSS couplings (for S = $ C ) in Table IVI 



cos 2/3. Note also that G and G+ C{V^H d H d ) and C(V^H U H U ) pieces 



f + is the usual pure left hand ed stau field with a mass 

™>f+ = VKJ 2 ~ m w 
are the neutral and charged Goldstone bosons that are 

absorbed by the Z and IF-bosons and are, therefore, the 

states with a zero eigenvalue in M Q and Mq, respectively. 

The physical states , $ e and <&o are related to the 

weak states $° E and $q y i a: 



derived from C^^L^L^) in (|56|) and from the following 



arc 



®E i = RE,ik^E,k , 



(65) 



OA 



Ro,ik^O,, 



where Rq, Re and Ro are the rotation matrices that 
diagonalize M c , M e and M Q , respectively. 

Notice that the mass matrices M c , M E and M Q de- 
pend only on four SUSY parameters: m A , m° au ,bz and 
tan j3. These parameters, therefore, completely fix the 
rotation matrices Rq, Ro and Re from which the CP- 
even and CP-odd neutral scalar spectrum as well as the 
charged scalar spectrum is completely determined. 

Clearly then, the type B RPV contributions involve 
the 3rd generation sleptons that can mix with the Higgs 
fields through a 63 bilinear RPV coupling which enters 
the slepton-Higgs mixed mass matrices in Here 
also, we can further sub-divide the type B RPV effects 
according to the type of scalar (S 1 ) and type of fermion 
(/) that are being exchanged in the loops: 



type Bl : S a = $ C , Q 
with a = 1, 3, i = 1, 2, 3 and 



ft 



type B2 : S a = $e, q and <F 



0,0 



fi 



(66) 



(67) 



C(V^H d H d ) = -i-H* d 



C(V,H U H U ) = -i-H* 



— r 3 ^ - —B, 

s w civ 



— r 3 Al + —B^ 
s w c w 



d%H d 



(68) 



(69) 

where the SU(2) scalar doublet fields L3, H d and H u are 
defined in (p0|). 



type Bl 



scalar (S a ) 
fermion 



2 j 

2 L(Zf) 
ij 

l R{Zf) 



°L{ a ) 



9z 



$c,a, a = 1, 3 
iti, i = 1,2, 3 

O-HZu) 

a-R(Zu) 



e m " a j ]^ a \^ 

V2s w m w s f3 C »J 



—e cot 28 w 8, 



Otf3 



TABLE VI: The couplings required for the calculation of 
T(Z — > didj) in the type Bl RPV scenario. The couplings 
follow the notation in (g), (0) and ([]). The couplings cll,r(Zu) 
are given in (H)-(H). 



with a = 1, 2, 3 , = 1, 3, i = 1, 2, 3. 

Note that we have omitted virtual exchanges of the 
two Goldstone bosons G and G + since the one-loop am- 
plitudes are being calculated in the unitary gauge. 

The two RPV effects, of types Bl and B2 above, are 
driven by the Higgs-slcpton scalar admixtures <&c, &e, 
$0 which couple to quarks through a combination of 
A' and Higgs Yukawa coupling. Hence, the Higgs-like 
components in <I>c, $e and «3>o will couple through the 
Higgs Yukawa terms, while the slepton-like component 
interact with the quarks via the A'-type RPV couplings 
in ((541). 



For the type B2 RPV case (see (|67|)) there are 10 one- 
loop diagrams that can potentially contribute to the de- 
cay Z — > djdj. These diagrams are depicted in Fig. [II]. 
The first eight diagrams in Fig. |l0| have the same topol- 
ogy as the generic diagrams of Fig. [I], while diagrams 9 
and 10 involve virtual exchanges of a Z-boson through 
the ZZ&e interaction. 

Using our generic notation for the one-loop amplitude 
in (j|), we calculate the form factors A I L J k , A^ k and 
Bjf k , -Bjj J fc with k = 1 — 10, which emerge from diagrams 
1-10 in Fig. [tOl (taking m di = 0, for 0?, s and 6-quarks): 
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<E>E 









FIG. 10: One-loop diagrams that contribute to the flavor changing decays Z — > didj in the type B2 RPV scenario. 



AlJ 

A L.\ 



AlJ 
A L,2 



AlJ 
A L,3 



AlJ 



AlJ 



A I. J 



AlJ 



9 a p h (E)iI h (0)iJ r l 



ot,{3,i 



9 B<x h (0)il ,{E)iJ r 2 



,{E)iI h {E)jJ 



aR(Zd) 2^ °L(a) 6 L(a) L 2C : 



-m z (C 23 — C 22 )] , 
-m| (C 23 — C| 2 )] , 

a L {Zd) E b( L(a) bi L(a) B l > 



h (0)il AO)iJ „7 
a L (Zd) 2^ °L(a) L(a) B l 



aV w = o 



(70) 
(71) 

(72) 

(73) 
(74) 

(75) 
(76) 



s £ J io — a -R(zd) 9zz^ E b L(d) [ 2 (^12 ~ Cn) 

a 

+J_(c' 10 + a i 1 1 )l . (79) 
m z \ 1 

Here also the right-handed form factors, A^ k and Bj^ k , 
are obtained from the corresponding left handed ones by 
interchanging L — > i? and i? — > L in all the couplings in 

The two-point and three-point loop form factors B\ 
with k = 5,7, C k x with a; G 11, 12, 21, 22, 23, 24 and k = 
1,2,3,4,9,10 and C* with a; e 0,11,12 and k = 9, 10 
which appear in (|7(i|)-(|79|) are given by: 



B\ = B 1 (m 2 d .,m% 0a ,m 2 dl S j , 



and 



(80) 
(81) 



where we have combined the contribution of the self 
energy diagrams 5+6 and 7+8: M.$ + Me = M^q 
and .M7 + .M8 = -Mrs, which leads accordingly to 

^L,5 + ^Lfi = ^L J 56 an d ^L J 7 + ^L^S = ^1/78- Also, 



B 



I J 



= for k = 1 - 



(77) 



nW _ 
B L,9 - 



a L(Zd)^2 9zZ<f E b L(J) J [ 2 (C12 - Cu) 
a 



(78) 



= Cx (m 2 d .,ml Ea ,ml O0 ,m dl ,q 2 ,m dj J , (82) 



£*2 q ( ^ 2 2 L1 2 



>™4,9 2 >™;L) , (83) 



C| = Ca, ^m| E a ,m^, m^,m^, g 2 ,^ ) , (84) 



c i = c x(ml 0a ,m 2 d .,m 2 d .,m 2 dj ,q 2 ,m 2 dl S j , (85) 
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2 2 2 







(86) 



C1U. y^*!U y^~* 
a: ' x — °i 



f 2 2 2 
Cr ^771^,777^, 777,^ 



2 2 2 



(87) 

where ^(mf , m^,p 2 ), C x {m(,mj,m%,p(,p%,p%) and 

C x {m\^m2,rn\ 1 p\ 1 p\,pl) are defined in the appendix. 

h (0)ij 

L(a) ' fi(a) ' U L(a) ' 

^R{a) ' Sz' 3 ano - .9zz$ E needed for evaluating the form fac- 



The couplings a L{zd) , a R{zd) , b 



tors above are given in Table VII. In particular, dL,R(Zd) 
are the SM left and right-handed couplings of the Z-boson 
to a pair of down quarks as given in (§-(§)■ The rest are 
obtained from the relevant interaction Lagrangian terms 
by rotating the SU(2) weak states &c,e,o to the phys- 
ical states <&c,e,o according to (pa). In particular, the 



§Edidj couplings b^^ a) and ^odid-j couplings b ( LR . t 



dO)ij 



follow the notation of the generic Sdf vertex in 
S = <& E or $o and f = d: 



i); for 



k{<S> E , a didj 



L(a) 



(89) 



These couplings emanate from both the Yukawa-like tri- 
linear RPV interactions in (^4|) and the neutral Higgs 
Yukawa vertices of a 2HDM of type II as given in (|23|). 

The coupling of a Z-boson to a $E,a®o,0 pair 
follow our generic definition of the VSS vertex i n (ffl ) . It 
is derived from the Lagrangian terms in ((5^), (pq) and 

The coupling gzz^ E °f ®E,a to a pair of Z-bosons is 
defined as: 



scalar (S a ) 
fermion (ft) 



~. z 3 

L(Zf) 
ij 

l R(Zf) 
h (E)ij 

u(E)ij 
°R(ac) 

iXO)ij 

°L( a ) 

iXO)ij 

°R( a ) 



al3 

9z 



9zz-s>, 



type B2 



&E, a , a = 1, 2, 3 and $o,a, a — 1,3 
di, i = 1,2,3 

a L(Zd) 
d-R(Zd) 



2s w mwC R E S ij I" — 2 ^i]i^E 



1 \l r>3a 



md i pln^ I 1 \'* p3a 



~ 1 2 SW m W c /3 R o' Si i + 72 A 3ji-RcT 
; e m rf^ pla; I ? \'* Z?3a 



sin 2» w \ u e n a — n E n Q + n E n Q 
-mz (cpR)? + Sf3 R 2 E a 



TABLE VII: The couplings required for the calculation of 
T{Z — > didj) in the type B2 RPV scenario. The couplings 
follow from the Feynman rules in © and @ 

The couplings dh,R(Zd) are given in (G 



Thus, without loss of generality, we trade the bilin- 
ear coupling 63 with a dimcnsionless RPV parame- 
ter, e, as follows: 



and is obtained from the following ZZ(® d and ZZ£® in- 
teraction terms (recall that £jj u are the CP-even SU(2) 
components of as defined in (|60|)) |22j |: 

£(^2,J = , 9g e , + ««e2) , (si) 

(2s w cw) 

where Vd and v u are the VEV's of the down-type and 
up-type Higgs doublets, respectively. 

Before presenting our numerical results for the type 
B RPV contribution let us discuss some of its salient 
features and outline the main assumptions and notations 
regarding the relevant parameter space involved: 

I. The pseudo-scalar "bare" mass (i.e., its mass in the 
RPC limit of b^ — > 0) can be approximated from 
the tree-level relation which, for t"i 3> 1, gives 
(m° A ) 2 ~ b tp 7 where b is the usual RPC soft- 
breaking bilinear Higgs term in the scalar potential, 
i.e., Vrpc D b H d H u . 



b 3 = e(m° A ) 2 cot [3 , (92) 

such that e ~ parametrizes the relative amount 
of RPV in the scalar potential. In particular, e <C 1 
for small bilinear RPV and e - 1 if RPV/RPC - 1 
in the SUSY scalar sector. 

The existing experimental limits on e come from: 

1. A non- vanishing 63 can generate a radiative 
(one-loop) tau-neutrino mass. The laboratory 
limit on the tau-neutrino mass allows, how- 
ever, the quantity ^ - e to be of - 0(1) @. 

2. The parameter 63, or equivalently the quan- 
tity e ~ can have important consequences 
on the CP-even and CP-odd Higgs-like scalar 
mass spectrum, see |36|, |37j]. In particular, 
e can drive the mass of the physical CP- 
even light Higgs below its present LEP2 lower 
bound which, for 777,^ ^ 200 GeV, is roughly 
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FIG. 11: BR(Z — > bs + bs) as a function of the "bare" sneutrino mass parameter m° sv (see text), for some combinations of 
values of rn\ and e (as indicated in the figure) and for tp = 3 (left plots) and = 50 (right plots). A332 = 1 is used and e is 
defined in (B2). 



m/j £ 110 GeV irrespective of tan/3 



[12] 



Also, a non-zero e can give rise to nega- 
tive eigenvalues (i.e., to the physical square 
masses) for the CP-even and CP-odd mass 
matrices Mf and Mq in and (||), de- 
pending on the values of the rest of the type 
B parameter space, i.e., on m A , mP sl/ and tp. 

Therefore, in what follows, we will vary the 
parameters {m A , m° sv , tp, e} subject to the ex- 
isting LEP2 lower bound on the light Higgs 
mass and to the requirement that mo^ , 
and m?+ are > 150 GeV. 

Since the light Higgs mass is very sensitive 
to higher order corrections to the 2x2 Higgs 
block in M§, as in |36|, |37| , in order to derive 
realistic exclusion regions for the parameter 
space {m A ,m° v ,tp,e} through the require- 
ment to/j ^110 GeV, we include the dominant 
higher order corrections (coming from the t — t 
sector) to the block in Mf., following 

the approximate formulae given in |39[l and 
taking the maximal mixing scenario with a 
typical squark mass of rriq = 1 TeV. 



[12] This bound is applicable in the maximal stop mixing scenario 
with a typical SUSY squark masses of 1 TeV. Note also that 
since i>3 / the hZZ coupling is smaller than its value in the 
RPC case leading to a smaller e + e _ — > Zh production rate. 
Thus, the limits reported in [ pq ] should be slightly relaxed in the 
type B RPV scenario. 



3. A non- vanishing e can also alter the cross- 
section for ZZ and WW pair production 
through s-channel exchanges of the CP-even 
scalars |36, |37|. The measured ZZ and 
WW cross-sections in LEP2 can thus be 
used to place limits on e as a function of 
{m A ,mP 8ll ,tp} . These limits, however, can be 
evaded if the v r ± e~ ~e + trilinear RPV coupling 
A131 is assumed small enough (see |3(], [37j). 
We will, therefore, not consider such limits be- 
low. 

II. Since 63 is not a flavor changing parameter, the 
transition between down-quarks of different gen- 
erations, i.e., between the external down-quarks 
di — > dj, is necessarily driven by a A' coupling 
with the appropriate non-diagonal indices (disre- 
garding flavor changing transitions due to small 
non-diagonal CKM elements). Thus, the type B 
RPV one-loop effect in Z — > didj is necessarily 
proportional to either b^X'^j j or b^X'^jj. 

In particular, for Z — > bs we find that the dominant 
contribution is attributed to the type Bl exchanges 
of the charged scalars and it arises when A 332 7^ 
0. The only other possible index combination for 
Z — > bs, which is A 323 ^ 0, yields a much smaller 
branching ratio. This enhancement for the (332) 
index combination can be traced to the fact that, 
for this particular combination, the charged scalar 
amplitude involves also a top-quark exchange, thus 
gaining a factor of mt /m c compared to the A 323 7^ 
case (which involves a charm-quark exchange in the 
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FIG. 12: Physical masses ol the heavy CP-even Higgs (rriff), CP-odd Higgs (mA), charged Higgs (m H +), CP-even tau-sneutrino 
(m+,), CP-odd tau-sneutrino [m~ v ) and the stau (m ST ), as a function of the "bare" tau-sneutrino mass (mjj„), for tp = 50, for 
m° A = 300, 600 or 900 GeV and for e = 0.4 (left figures) and e = 0.9 (right figures), e is defined in (§§). 



loops). 

III. In the limit e — > the type B2 effect vanishes. How- 

ever, since e — > causes the charged Higgs sector 
to decouple from the stau sector and since the RPC 
MSSM Higgs sector is similar to the 2HDM of type 
II, the type Bl contribution approaches that of the 
type II 2HDM in this limit. Thus, for e -> 0, the 
type Bl RPV effect will be proportional to the off- 
diagonal CKM elements as in the case of the type 
II 2HDM discussed in section III. 

IV. In the numerical analysis below we will set A 332 = 

1, while all other lambda's with different index 
combinations are set to zero. The experimental 
limit on this coupling, derived from i?< > = T(Z — > 
hadrons)/r(Z -> ££) M, is (at the 2a level) 



:s:;2 



0.45 for squark masses of ~ 100 GeV , 
while A 332 = 1 is allowed for squark masses £ 650 
GeV. The perturbativity bound on this coupling 
is A 332 = 1.04 Thus, we will assume that 

the squarks are heavy enough to allow A 332 to lie 
near its perturbativity limit (recall that no squarks 
are involved in the type B RPV contribution to 
Z^bs). 

In Fig. [ll] we show BR(Z — ► bs + bs) as a function of 
the "bare" tau-sneutrino mass vnP sv (i.e., what would be 
its mass in the RPC limit), for various possible values of 
mP A , e and for tp = 3 (left side) and tp = 50 (right side) 



p4j . Evidently, BR(Z — » bs + bs) is much larger in the 
high tan f3 scenario and it drops with m° sl/ . 

The masses of the heavy CP-even Higgs, CP-odd Higgs 
and charged Higgs as well as the CP-even, CP-odd tau- 
sneutrino and the stau particles are depicted in Fig. |lj, 
for tp = 50 and for the same combinations of s and 

We note that in the 



that are used in Fig. 
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limit (jtia) ^ m z (applicable to the values of m A in 
Figs. O and |l2| ) one has run ~ ~ mH+ and if 
in addition (m^) 2 3> m%, then also the CP-even, CP- 
odd tau-sneutrinos and the stau are roughly degenerate. 
Thus, only two curves are shown in each plot in Figs, [l^, 
which arc sufficient to approximately describe all these 
six scalar masses. 



Fig. 12 shows that at some instances, the Higgs-like 
and slepton-like scalar masses exhibit a discontinuous 
jump, at which point they "switch" identities. This phe- 
nomena is caused by the particular dependence of the 
physical scalar masses on the "bare" masses m A and 
m°„ in the presence of e ^ 0. In particular, the correc- 
tions to the "bare" scalar masses due to a non- vanishing 
63 term are proportional to factors of [(m^) — (m^)] 
(for more details see [|6[ |37j]), thereby changing sign at 
the turning points. Moreover, the off-diagonal elements 
of the rotation matrices Re, Ro and Rc-i which are re- 
sponsible for the slepton-Higgs mixings, are also inversely 
proportional to factors of [(W4) — (^L)] > therefore, en- 
hancing the type B RPV effect as m A approaches m° sv as 
can be seen in Figs. [ll]. 
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model 


scalars in the loops 


BR(Z -> 6s) 


SM 


VF-boson (no scalars) 


1CT 8 


2HDMII 


charged Higgs 


kt 10 


T2HDM 


charged Higgs 


KT 8 


SUSY with t — c mixing 


t — 5 admixtures 


KT 8 


SUSY with b — s mixing 


6 — s admixtures 


KT 6 


SUSY with trilinear J?p-violation 


squarks and sleptons 


10 -io 


SUSY with bilinear i?p-violation 


slepton-Higgs admixtures 


KT 6 



TABLE VIII: The best case values for the branching ratio of Z — + 6s for each of the six models considered in this paper upon 
imposing the available experimental limits on the relevant parameter space of each of them. The SM prediction is also given. 



To summarize this section, with a large tan/3, a 
BR(Z -> 6s + bs) ~ O(10" 6 ) is possible within the 
type B scenario, e.g., for 40% lepton number violation in 
the SUSY scalar potential (e = 0.4) and if the sleptons 
masses lie around ~ 200 GeV. For a heavier slepton spec- 
trum a larger e is required in order to push the branching 
ratio to the 10 -6 level. 

It should also be emphasized that since BR(Z bs + 
bs) is dominated by the A 332 , the decay Z — > bs is an 
efficient probe of this specific flavor changing trilinear 
RPV coupling. 



VI. EXPERIMENTAL FEASIBILITY 

In this section we will very briefly comment about the 
feasibility of observing (or achieving a limit) a signal of 
Z — > bs with a branching ratio of order 10" 6 , at a Linear 
Collider producing 10 9 Z-bosons. 

Such a signal should appear in the detector as an event 
with one b-jet and one light-jet (assuming no distinction 
is made between light — d, u or s-quarks). In the spirit 
of the analysis made with the 1993 and 1994 LEP data 
on Z — > bs ||, one defines and to be the efficien- 
cies that a quark (or anti-quark) of flavor q is tagged 
as a 6- jet (B) and light-jet (L), respectively. Thus, the 
key efficiency parameters for the detection of Z — > bs are 



fied as Z — > B + L, is: 



c fe ' Might 



and e^, where the latter represent the proba- 
bility that a 6-jet is identified as a light-jet and is im- 
portant for controlling the dominant background (to the 
Z — ► B + L signal) coming from Z — > 66. Note, that 
due to the expected smallness of the " Purity" parameters 
£jf ght , and (see ||), the background to Z — > B + L 
caused by the SM Z — > dd, ss, uu, cc decays will be 
sub-dominant. 

With 10 9 Z-bosons, the expected number of events 
coming from Z — > bs (i.e., from new physics) and identi- 



S - 10 9 x e£eL u BR(Z -> 6s) . 



(93) 



Similarly, the expected number of background Z — * B+L 
events coming from the SM decay Z — > 66 is: 



B - 10 J 



f B f L 



BR(Z -> 66) 



(94) 



Using (|93|) and (|94[), the expected statistical signifi- 
cance, S/vB, of the new physics signal Z — * 6s, with 
a branching ratio of order 10 -6 , can reach beyond the 



3-sigma level for eir ~ 0.6 — 0.8, e 



light 



0.3 - 0.5 and 



~ O(10 4 ). These values require an improvement to 
the 1993 and 1994 analysis [§, by a factor of 2-3 for ef 
and e^ght and by an order of magnitude for e^. With 
the expected advancement in the jet-tagging methods, in 
particular, for two-body decays of the Z-boson, these re- 
quired values for the efficiency parameters above should 
be well within the reach of the future Linear Collider. 

We can also get a clue about how low one can go in the 
value (or limit) of BR(Z — > bs) with 10 9 Z-bosons, from 
the fact that the LEP preliminary results || achieved 
BR(Z -> 6s) < O(10- 3 ) with O(10 6 ) Z-bosons. Scal- 
ing this limit, especially with the expected advance in 
b-tagging and identification of non-b jets methods, an 
O(10 -6 ) branching ratios should be easily attained at a 
Giga-Z factory. 



VII. SUMMARY 

We have re-examined the flavor changing radiative de- 
cays of a Z-boson to a pair of down-quarks, Z — > didj, 
with I J. These Z-decay channels may prove useful 
in searching for new flavor physics beyond the SM at the 



25 



TESLA collider, or any other future collider, which may 
be designed to run on the Z-pole with high luminosi- 
ties, thus accumulating more than 10 9 on-shell Z-bosons. 
With advances in technology, e.g., improved b-tagging ef- 
ficiencies, the flavor changing decay Z —> bs - most likely 
the easiest to detect among the flavor changing hadronic 
Z-decays - may be accessible to a Giga-Z option even for 
branching ratios as small as BR(Z — > 6s) ~ 10~ 7 — 10~ 6 . 

The di — > dj transition was assumed to be generated 
at one-loop through flavor violation in interactions be- 
tween scalars and fermions. 

A complete analytical derivation of the width T(Z — > 
didj) is presented using the form factor approach for the 
Zd]dj interaction vertex. These form factors are evalu- 
ated for the complete set of scalar-fermion one-loop ex- 
changes with generic scalar-fermion flavor-violating cou- 
plings. 

This prescription is then applied to the decay Z — > bs 
in six beyond the SM model scenarios for flavor-violation 
in the scalar sector: 

1. Two Higgs doublet models with non-standard 
charged-Higgs couplings to quarks: 

• A two Higgs doublet model of type II (2HD- 
MII). 

• A two Higgs doublet model "for the top- 
quark" (T2HDM). 

2. Supersymmetry with flavor- violation in the squark 
sector: 

• Supersymmetry with stop-scharm mixing. 

• Supersymmetry with sbottom-sstrange mix- 
ing. 

3. Supersymmetry with flavor-violation from R-parity 
violating interactions: 

• Supersymmetry with trilinear R-parity viola- 
tion. 

• Supersymmetry with trilinear and bilinear R- 
parity violation. 



Folding in the existing experimental limits on the rel- 
evant parameter space of each of these models, we cal- 
culated the branching ratio for the decay Z — > bs. The 
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highlights of our results are summarized in Table 
In particular, we find that two Higgs doublet models with 
flavor violation originating from charged scalar interac- 
tions with fermions are expected to yield an extremely 
small BR(Z — > bs); smaller than the SM prediction and 
smaller than the reach of a Giga-Z £ + £~ collider. Thus, a 
signal of Z — > bs in TESLA will be inconsistent with the 
underlying mechanisms for flavor violation in these two 
Higgs doublets model and will, therefore, rule out these 
options. 

The same conclusions can be drawn in the stop-scharm 
mixing and the trilinear R-parity violation SUSY scenar- 
ios. On the other hand, SUSY with mixings between 
the bottom and strange- type squarks and/or mixings be- 
tween sleptons and Higgs fields (bilinear R-parity viola- 
tion) both of which may originate from the soft SUSY 
breaking sector, can drive the BR(Z — » bs) to the 10~ 6 
level for large tan (3 values. This enhancement is typi- 
cal to these two flavor-violating SUSY scenarios if there 
are large mass-splittings between the scalars exchanged 
in the loops due to a GIM-likc cancellation which is oper- 
ational in the scalar mass-matrices and is, therefore, less 
effective as the scalar masses depart from degeneracy. 

A Z — > bs signal in a Giga-Z TESLA or any other col- 
lider may, therefore, be a good indication for the under- 
lying dynamics of these two flavor-violating SUSY sce- 
narios and, if interpreted in that way, will provide for 
evidence of an hierarchical structure in the mass spec- 
trum of the SUSY scalar sector. 
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APPENDIX A: ONE-LOOP FORM FACTORS 

In this appendix we give the two-point and three-point 
one-loop form factors which are defined by the one-loop 
momentum integrals as follows E3] : 



J 



C ; C M ; C^ u (ml,ml,mj,pl,pl,pf) 



d 4 q 



in 2 [q 2 - m\] [(q + pi) 2 ~ mf] [(q - p 3 ) 5 



(Al) 
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fi ri ( 2 2 2 2 2 2\ 

6 ; C p (m 1 ,m 2! m 3 ,p 1 ,p 2 ,p 3 ) 
B Q ; Bfj, (ml,ml,p 2 ) 



d 4 q 



g 2 ; g 2 g M 



27T 2 [q 2 



m i] [(« + pi) 2 - [(<? - P3) 2 - 
i; <?m 



Z7T 2 [g 2 — to 2 ] [(q +p) 2 — m|] 

I 



where '^2 i Pi — is to be understood above. 

The coefficients B x with x G 0,1, with x S 
0, 11, 12, 21, 22, 23, 24 and C x with a; £ 0, 11, 12 are then 



defined through the following relations 43 



Bf, = p^Bt , 



(A4) 
(A5) 



Plti.Cn + P2tiCl 



and 



(A2) 
(A3) 



(A6) 



C^ v =Pl l _ l PlyC 2 l +P2^P2vC 2 2 + {VlVl}^ C 2 3 +9^C 2 4 , 

(A7) 

where {ab}^ = a^b v + a u b^. 



